Biochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms

The actin polymerization machinery, comprising the ARP2/3 complex and its activators, the WASP family proteins, has been implicated in regulating a broad spectrum of nuclear processes, such as transcriptional regulation and nuclear organization. Here, using clustered protocadherin (cPcdh) and {beta}-globin genes as model systems, we showed that WAVE2, a member of the WASP family, regulates chromatin organization by maintaining heterochromatin dynamics. Specifically, by CRISPR DNA-fragment editing, in conjunction with integrated analyses of ChIP-seq, MeDIP-seq, ATAC-seq, 4C-seq, and RNA-seq, we showed that deposition of H3K9me3, a key heterochromatin mark, is significantly decreased at the cPcdh locus upon WAVE2 deletion, concurrent with aberrant accumulation of CTCF/cohesin complex at promoter regions and spatial reorganization of chromatin architecture around nucleolus. In addition, REST/NRSF exerts a similar heterochromatindependent effect on the cPcdh locus. Finally, genetic and genomic data showed that WAVE2 regulates {beta}-globin gene expression by maintaining heterochromatin status. Together our data suggested that WAVE2 and REST/NRSF regulate clustered gene expression in a heterochromatin-dependent manner.

The cockroach Blattella germanica possesses panoistic ovaries, in which oocytes lack nurse cells and therefore need to rely on their own transcriptional activity to support embryogenesis. Ovarian development in this species involves the development of a single basal ovarian follicle (BOF) per gonadotropic cycle, a process strictly regulated by endocrine signals, primarily juvenile hormone and ecdysone, which act at both the transcriptional and translational levels. In addition, transcriptional activity in these ovaries is necessary for both regulating and genome protection, and at this level, PIWI-interacting RNAs (piRNAs) play an essential role. Although insect ovaries are known to be particularly rich in piRNAs, their function in ovary maturation is still not well defined. For this purpose, we characterize the piRNA expression dynamics across seven key developmental and reproductive stages, ranging from late nymphal instars to post-vitellogenic adults. piRNA expression in B. germanica shows coordinated fluctuations. Expression remains stable in previtellogenic ovaries, whereas vitellogenic ovaries show pronounced changes. Moreover, vitellogenic ovaries exhibit reduced piRNA diversity due to strong enrichment of a subset of highly expressed piRNAs. Our data show that although piRNAs predominantly map to transposable elements, particularly LINEs, there is a notable increase in gene-derived piRNAs toward the end of the cycle. Our results suggest regulatory roles of piRNAs in modulating both TEs and mRNAs during BOF maturation, likely related to changes in the follicular cell program.

Muscleblind-like (MBNL) RNA-binding proteins (RBPs) possess modular domains that mediate regulation of alternative splicing and RNA localization. Myotonic Dystrophy Type 1 is a CTG repeat expansion disorder where MBNL is sequestered into intranuclear RNA foci, impairing its function. Previous studies found that MBNL self-associates through its exon 7, but the nature of this interaction is not well understood. We identified a cysteine in MBNL1 exon 7 that enables dimerization through formation of an intermolecular disulfide bond. We likewise demonstrate that MBNL2 dimerizes by forming disulfide bonds between multiple cysteines in its carboxy-terminus. Nucleocytoplasmic fractionation revealed a greater proportion of MBNL1 dimer in the nucleus, suggesting a nuclear function for the MBNL1 dimer. We investigated a connection between MBNL1 dimerization and MBNL1-mediated regulation of alternative splicing. To accomplish this, we mutated the MBNL1 cysteine in question to alanine (C325A) and performed RNAseq. We uncovered novel splicing events sensitive to MBNL1 dimerization. We also found that MBNL1 C325A, when co-expressed with expanded CTG repeats, produces smaller, more numerous foci, suggesting a role for the MBNL1 dimer in maintaining foci integrity. These results provide insight into biological and pathological mechanisms of MBNL1 dimerization and suggest other RBPs might similarly dimerize to regulate function. GRAPHICAL ABSTRACT

Enhancer RNAs (eRNAs) are non-coding transcripts produced at enhancer regions, which appear to be involved in transcriptional regulation. Up to date, these have been primarily investigated using labor-and cost-intensive genomic techniques. However, the precise mechanisms by which eRNA transcription or the eRNA transcripts themselves mediate transcriptional regulation remain unclear. Here, we present a novel experimental approach that allows us to analyze the characteristics of eRNA transcription in fixed and live whole Drosophila melanogaster embryos. We employ the anterior-posterior patterning genes as a model system to investigate the dynamics of eRNA expression, utilizing an imaging-based approach. We combined high-sensitivity fluorescence in situ hybridization (FISH) chain reaction (HCR) with high-resolution confocal microscopy to detect eRNA and mRNA molecules. Through this experimental assay, we identified foci of elevated transcriptional activity that generate eRNA transcripts correlated with mRNA production at the same gene locus. We could show that this eRNA transcription is independent of promoter activity. Additionally, we demonstrate that insulators can influence eRNA transcription, resulting in loss of eRNA transcription. Moreover, we observe that eRNAs can originate both within classical enhancer regions and outside of them, including from foreign bacterial sequences when these are placed near enhancer sequences, underscoring the strong influence of local regulatory context on eRNA initiation. In live embryos using MS2-MCP live imaging, our analysis of insulators showed a modest reduction in mRNA burst intensity accompanied by a slight increase in burst frequency. Overall, our imaging-based approach offers a novel platform for dissecting enhancer-eRNA interactions and could be adapted for wider applications.

The ubiquitin conjugating enzyme UBE2J1/Ubc6e localizes to the endoplasmic reticulum where it mediates the ubiquitination and proteasomal degradation of terminally misfolded proteins. Although the protein is known to undergo phosphorylation at serine S184, we have considered modification at an additional site and used a bespoke anti-phospho antibody to confirm phosphorylation also at serine residue S266. Despite the well-described role of UBE2J1 in ER associated degradation (ERAD), we found no evidence for regulation at S266 during Unfolded Protein Response (UPR) induction by thapsigargin. Instead, our studies suggest that phosphorylation occurs independently at the S184 and S266 sites, with mutation at one site failing to disrupt basal phosphorylation at the second. We identified several contexts in which these two phosphorylations were differentially regulated. For example, ER localization, which is important for phosphorylation at S184, was not required for modification at S266, and sensitivity to proteasome inhibitors, which is regarded as a distinguishing feature of the S184 phospho-variant, was unaltered by the S266A mutation. Regarding regulation at S266 on the other hand, we found that pharmacological activation of protein kinase A resulted in rapid phosphorylation, with differential use of phospho-specific antibodies confirming that phosphorylation at S184 was unchanged by this treatment. Hormonal stimulation by glucagon resulted in a similar pattern of UBE2J1 phosphorylation, which occurred exclusively at S266 and could be inhibited by H89. The differential regulation demonstrated in these studies extends our understanding of the UBE2J1 enzyme, and may indicate a role in the integration of energy metabolism with environmental stress conditions.

Germ cells are the only cells of an organism that pass onto the next generation and, hence perpetuate the species. To ensure this, germ cells need dedicated transcriptional repertoire, that ensure specification, proliferation, differentiation and fate maintenance. We previously characterized LSL-1, a conserved zinc-finger transcription factor that acts as a major direct transcriptional activator of genes involved in germ cell development, fate specification, meiosis and genome stability. Here, we show that LSL-1 interacts with the transcription factor HIM-17, the chromatin proteins BRA-2 and XND-1. These proteins are functionally related to LSL-1 and they colocalize at germline gene promoters, forming most likely a transcription-promoting complex. Furthermore, LSL-1 lies in close proximity to members of the COMPASS and the MOF complexes, corroborating the observation that HIM-17 and LSL-1 are required to maintain normal level of H3K4 methylation in the gonad. Finally, we show that LSL-1 interacting partners are necessary to maintain germ cell fate. Altogether, we propose that LSL-1 interacts with transcription regulators and chromatin modifiers to ensure the establishment of the transcriptional repertoire appropriate for germline function as well as for cell fate maintenance.

MicroRNAs (miRNAs) play essential roles in the posttranscriptional regulation of gene expression in organisms. In the process of synthesizing mature miRNAs from miRNA precursors, the miRNA precursors are cleaved via Dicer at their loop structure, after which the miRNA precursors become mature and regulate transcription. However, the consequences of altering the loop sequence are not fully understood. The silkworm Bombyx mori is a lepidopteran insect with many genetic strains. We identified a mutant of the miRNA miR-3260 whose the part of the loop structure was lacking in a silkworm strain with translucent larval skin. Here, we aimed to analyze the role of wild-type miR-3260 and the influence of the mutation of the loop structure in B. mori. First, we identified the genomic region responsible for the translucent larval skin phenotype and determined that the mutated miR-3260 nucleotide sequences. Then, we predicted the binding partners of wild-type miR-3260 using the RNA hybrid tool and found two juvenile hormone (JH)-related genes as targets of wild-type miR-3260. Next, we assessed the relationships between miR-3260 and JH and found that miR-3260 was highly expressed in the Corpora allata and its expression responded to JH treatment. Meanwhile, miR-3260 mimic and inhibitor did not induce the typical phenotypes associated with JH in B. mori. Then, we compared the dicing products from wild-type and mutant miR-3260 precursors and observed that neither form underwent Dicer-mediated cleavage when the loop structure was altered. These results suggest that loop mutations in the miR-3260 precursor may not influence dicing activity, consistent with the lack of observable phenotypic effects.

Caffeine, the most widely consumed stimulant worldwide and primarily sourced from coffee, is well known for its central nervous system effects. Emerging evidence indicates that caffeine also modulates key cellular processes, including DNA repair. It inhibits the kinase activity of ATM and ATR-essential DNA damage response proteins, and impairs homologous recombination (HR)-mediated repair through multiple mechanisms. However, its effects on nonhomologous end joining (NHEJ), a major double-strand break (DSB) repair pathway, have been underexplored. In a recent study, we reported that caffeine inhibits NHEJ primarily by interfering with Ligase IV/XRCC4 complex, using in vitro and ex vivo model systems. Given coffees role as a primary dietary caffeine source, this study investigates the impact of Coffea arabica decoction on NHEJ-mediated DSB repair. High-performance liquid chromatography (HPLC) quantified caffeine levels in the decoction, followed by in vitro and ex vivo assays to evaluate NHEJ efficiency. Results demonstrate that coffee decoction inhibits end joining of both compatible and noncompatible DNA ends in cell-free systems derived from normal and cancer cells. Extrachromosomal repair assays confirmed impaired intracellular NHEJ, leading to accumulation of unrepaired DSBs in human cells. Kinetic analysis of {gamma}-H2AX foci formation and resolution revealed persistent DNA breaks and reduced repair kinetics. Reconstitution experiments verified that the decoction specifically targets the Ligase IV/XRCC4 complex. These findings, building on our previous work, establish coffee decoction as a potent NHEJ inhibitor, mirroring purified caffeines effects. This underscores caffeines interference with endogenous DNA repair, with profound implications for cancer therapy by sensitizing tumors to genotoxic treatments.

Heavy metal ATPases (HMAs) are important group of transmembrane proteins involved in homeostasis of metal ions in plant systems. In this study, a comprehensive analysis of genome assembly (VC1973A v7.1) resulted in the identification of nine HMA genes (VrHMA) and their corresponding proteins in Mungbean, an agronomically important legume crop known for its nutritional values. VrHMA proteins were also characterized based on their biomolecular features, conserved domains and motifs arrangement, transmembrane helices, pore-line helices, subcellular location and occurrence of signal peptides. Based on sequence homology, nine VrHMAs were clustered into two major substrate-specific groups: VrHMA1, VrHMA5 and VrHMA7 were categorized under the Zn/Co/Cd/Pb ATPase group, whereas the remaining six VrHMAs belong to the Cu/Ag subgroup. Gene structure analysis and promoter scanning revealed the structural divergence and presence of various stress-responsive cis-acting elements, respectively. The expression analysis of VrHMA genes in root and leaf tissues, in response to heavy metal (Zn, Cd and Cu) stress, indicates their role in the uptake, transport and sequestration of metal ions. Interestingly, VrHMA5 showed incremental upregulation in roots in response to all three heavy metal stresses, whereas its expression was only upregulated in the leaf tissues under Zn stress, which indicates its role in vascular transport in V. radiata. In addition, this study provides valuable insights into the functional roles of VrHMA genes and will lay a foundation for future genetic improvement in mung bean aimed at enhanced heavy metal stress tolerance and micronutrient homeostasis.

BackgroundPreeclampsia, a maternal hypertensive syndrome affect fetal brain development and cerebral angiogenesis, with potential acute and long-term consequences. Underlying mechanisms of these brain vascular alterations are unknown. This study investigates the role of thrombospondin-1 (TSP-1), an antiangiogenic glycoprotein, as a key mediator of communication between the fetoplacental and fetal brain endothelium in the context of preeclampsia. MethodsConditioned media (CM) of human umbilical vein endothelial cells (HUVECs) from normal pregnancies (NP-CM) and preeclamptic pregnancies (PE-CM), were used to treat human (hCMEC/D3) and murine brain microvascular endothelial cells (BMECs). A proteomic analysis was performed in plasma of the umbilical cord of normal pregnancy and preeclampsia. TSP-1 was identify using proteomic analysis and confirmed by Western blot. PE-CM depleted of TSP-1, using immunoprecipitation, was used to evaluate protein-protein interaction with vascular endothelial growth factor (VEGF). Antibody-mediated blockage of TSP-1 was used to investigate antiangiogenic effect and pro-angiogenic signaling pathways in brain endothelial cells exposed to PE-CM. ResultsPE-CM significantly reduced angiogenesis, migration, and invasion of brain endothelial cells and altered cytoskeletal organization. These effects were accompanied by reduced VEGFR2 and AKT signaling, indicating impaired angiogenic pathways. Proteomic analysis of umbilical cord plasma revealed elevated TSP-1 levels in preeclampsia, which was confirmed by Western blotting. TSP-1 was also increased in PE-CM, and immunoprecipitation assays suggested a protein-protein interaction with VEGF. Antibody-mediated blockade of TSP-1 restored angiogenesis, as reflected by increased total tube length, and rescued VEGFR2 and AKT signaling in brain endothelial cells exposed to PE-CM. ConclusionTSP-1-mediated endothelium-endothelium communication between placenta-brain axis in offspring of mothers with preeclampsia. This communication mediated by TSP-1 may contribute to acute and long-lasting cerebrovascular dysfunction observed in infants exposed to preeclampsia.

Astrocytes play a key role in neuronal homeostasis and in various neural disorders. The generation of astrocytes from neural progenitor cells (NPCs) and its functions are under a complex control of several signaling networks and transcription factors. In this study, we demonstrate that the transcription factor, GLIS similar 3 (GLIS3), which has been implicated in several neurodegenerative diseases, is highly expressed in astrocytes, and is required for the efficient differentiation of human NPCs into astrocytes. Loss of GLIS3 function greatly impairs astrocytes differentiation, resulting in reduced expression of astrocyte markers, whereas expression of exogenous GLIS3 restores the induction of astrocyte specific genes indicating a critical role for GLIS3 in astrocyte differentiation. Integrated transcriptomic and cistromic analyses revealed that GLIS3 directly regulates the transcription of several astrocyte-associated genes, including GFAP, SLC1A2, NFIA, and ATF3, in coordination with lineage-determining factors, such as STAT3, NFIA, and SOX9. We hypothesize that GLIS3 dysfunction disrupts this transcriptional network thereby contributing to astrocyte-associated neurological disorders. Identification of GLIS3 as a key regulator of astrocyte differentiation and gene expression will advance our understanding of its role in neurodegenerative diseases and may provide a new therapeutic target.

Osteosarcoma (OS) is an aggressive bone cancer that most commonly affects children and young adults. OS exhibits a high degree of genomic complexity, as well as cellular plasticity, and dynamic transcriptional regulation is suggested to contribute to treatment resistance and metastasis. Cell lines are well characterized as models to advance our knowledge on OS biology. HOS and U2OS cells have increased invasiveness and higher migratory ability compared with MG63. In this study, we employed a tandem array of consensus transcription factor response elements (catTFREs) proteomic approach to characterize transcription factor (TF) regulatory networks related to OS aggressiveness. We mapped 7,594 proteins and enriched 352 transcription factors and coregulators. When we integrated proteomics with cell line specific gene expression and chromatin accessibility we classified the proteins into different OS cell line dependent sub-clusters and identified TFs and coregulators common for all cell lines and specific for individual cell lines. We demonstrate that RUNX2, MYBL2 and HMGA2 are specifically enriched in HOS and U2OS and may be linked to the cell aggressiveness. ETV5, JUNB, NFIX and ZEB1 were among TFs specific to MG63. Our analysis provides a more comprehensive understanding of the transcriptional drivers that shape OS regulatory landscapes and may have future therapeutic implications.

Long noncoding RNAs (LncRNAs) have emerged as a class of important regulatory ncRNAs and are known to fine-tune numerous cellular processes including proliferation, differentiation and development; however, their role in quiescence still remains largely unexplored. A miRNA host gene lncRNA, MIR503HG, has been reported to play important role in cancer development. Here, we demonstrate the role of MIR503HG lncRNA in regulating cellular quiescence. MIR503HG displays elevated levels in human diploid fibroblasts induced to undergo quiescence. Depletion of MIR503HG in HDFs affects the entry of cells into quiescence but has no effect on cell cycle progression, suggesting its role in quiescence attainment and/or maintenance. Additionally, MIR503HG depletion led to a drastic decrease in the levels of miR508 target, PTEN with a concomitant increase in pAkt levels, indicating its role in negative regulation of miR508. Further, we demonstrate that the lncRNA MIR503HG regulates PTEN levels by acting as a ceRNA for miR508 to maintain cellular quiescence. Our studies illustrate that MIR503HG can function synergistically with miR503 to maintain cells under quiescence and both the miRNA-HG and the miRNA encoded by its gene locus synergistically control the same biological process in different ways by regulating different downstream genes.

Pre-mRNA splicing is an essential step in eukaryotic gene expression during which spliceosomes remove introns from nascent RNAs while ligating the adjacent exons. Spliceosomes are cellular nanomachines composed of five small nuclear (snRNA) components and dozens of proteins, most of which are highly conserved. Despite the high conservation of many splicing factors between S. cerevisiae and H. sapiens, several protein components of the S. cerevisiae spliceosome are not essential for growth under normal laboratory conditions. This is particularly surprising for nonessential factors whose conserved domains contact the spliceosomes catalytic core. Uncovering a function for these splicing factors can be challenging since they are not required for viability, may engage in functionally redundant interactions, and may display only weak phenotypes in the absence of secondary mutations in other spliceosome components. One such nonessential factor is the Cwc15 protein. Cwc15s highly conserved N-terminus directly contacts the U2/U6 di-snRNA within the spliceosome catalytic core; yet its precise role in splicing has not been defined in any organism. In this work, we use molecular genetics in S. cerevisiae combined with splicing reporter assays to study Cwc15p function. We propose that Cwc15p not only promotes active site stability during 5 splice site cleavage but also impacts structural transitions into and out of this spliceosome conformation. This function may be critical for splicing in S. cerevisiae under nonoptimal conditions, facilitating use of weak or alternate splice sites, and could have implications for proofreading of spliceosome active site formation. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=146 SRC="FIGDIR/small/713263v1_ufig1.gif" ALT="Figure 1"> View larger version (74K): org.highwire.dtl.DTLVardef@b296c5org.highwire.dtl.DTLVardef@c87b91org.highwire.dtl.DTLVardef@287011org.highwire.dtl.DTLVardef@d59741_HPS_FORMAT_FIGEXP M_FIG C_FIG Article SummaryPre-mRNA splicing is carried out by large macromolecular machines called spliceosomes which are composed of several snRNAs and dozens of proteins. Despite decades of study, the functions of many splicing factors such as S. cerevisiae Cwc15p remain unknown. Cwc15p is highly conserved among eukaryotes and directly contacts the spliceosome catalytic core. Here, we have used genetic and splicing reporter assays to study the function of Cwc15p during splicing in vivo. We propose that Cwc15p both stabilizes the spliceosome active site during 5 splice site cleavage and impacts remodeling of that site.

Mammalian telomeric DNA comprises long tracts of tandem TTAGGG repeats. The same repeats are also found at internal chromosomal regions called interstitial telomeric sequences (ITSs). Telomeres are transcribed into UUAGGG-containing transcripts, named TERRA, which serve multiple functions in maintaining telomere integrity. Complementary RNAs containing C-rich telomeric repeats, named ARIA, have also been identified in few yeast mutants and mammalian cells with dysfunctional telomeres. The molecular features and functions of ARIA remain understudied, mainly due to its low abundance and the lack of suitable cellular systems. Here, we show that Chinese hamster ovary (CHO) cells produce abundant TERRA and ARIA transcripts, predominantly originating from ITSs. Both RNAs are polyadenylated, exhibit relatively short half-lives and form large cellular foci. We also show that ARIA depletion leads to exposure of single-stranded (ss) DNA at ITSs and that ssDNA exposure increases when ITS DNA is damaged. SsDNA formation does not require the DNA damage signaling kinases ATM and ATR, nor the exonucleases DNA2 and EXO1; however, ATM prevents excessive ssDNA accumulation when ARIA function is inhibited. These findings establish CHO cells as a powerful model to dissect telomeric RNA functions and reveal ARIA as a key regulator of telomeric repeat DNA integrity.

Ehmt2 is a key H3K9 methyltransferase that regulates genome silencing and structural integrity during animal development. In addition to this canonical function, Ehmt2 has also been implicated in neural tissues mediating both direct and indirect transcriptional activation, and exon splicing, to facilitate proper neural cell differentiation and survival. Several germline loss-of-function animal models have been developed showing both conserved and divergent phenotypes that range from embryonic lethality to behavioural deficits in adult, fertile animals. Here, we generated the first maternal-zygotic ehmt2 loss of function mutant in zebrafish using CRISPR-Cas9 mutagenesis. An assessment of the pattern of H3K9 methylation in mutant embryos by ChIP-seq indicates that there are aberrant levels of this repressive mark, including reduction in discrete 5 non-coding regions of genes, but with no significant change in the overall pattern distribution of these marks across the genome. Global transcriptome and morphological analyses demonstrated that mutant embryos displayed greater variation in the timing of developmental progression that is, on average, slower compared to controls. Despite this, mutant embryos ultimately survive and are fertile. Through examination of progenitor cell dynamics and gene expression profiles, we found that the delay in embryonic development was associated with longer rates of S-M phases of the progenitor cell cycle in mutants leading to deficits in tissue growth. Finally, our data suggest a robust network of epigenetic regulators can potentially compensate for Ehmt2 loss of function and permit embryonic development and survival in ehmt2 mutant zebrafish. Our work establishes a zebrafish ehmt2 loss of function model that will facilitate examination of the complex and varied roles of Ehmt2 in vertebrate development.

Mechanosensing involves proteins detecting mechanical changes in the cytoskeleton or at cell adhesion sites. These interactions initiate signaling cascades that produce biochemical effects such as post-translational modifications or cytoskeletal rearrangements. Filamin is a ubiquitous mechanosensing protein that binds actin filaments and senses pulling forces within the cytoskeleton. Drosophila Filamin (Cheerio) is structurally similar to mammalian Filamin, with roles in egg chamber development, embryo cellularization, and integrity of muscle attachment sites and Z discs in Drosophila indirect flight muscles (IFMs). Here we report a potential novel binding partner of Drosophila Filamins: the death-associated protein kinase Drak that functions as a myosin light chain kinase. We found that Drak biochemically bound to an open mutant of Filamin that resembles the mechanically activated form partially bound to wild type Filamin and did not bind to closed mutant of Filamin. The interaction site was mapped to the intrinsically unfolded C-terminal region of Drak. To study the functional role of Drak-Filamin interaction, we studied two developmental events where Drak has been earlier shown to be expressed and where Filamin also functions: early embryonic cellularization and indirect flight muscle development at pupal stages. We found partial colocalization between Drak-GFP and Filamin-mCherry during the initiation of cellularization furrow, and at the time of myotube attachment site maturation in tendon cells. However, functionally we could not show direct correlation between Filamin and Drak. Our studies reveal interesting new expression patterns of Drak during Drosophila development and provide detailed information about Filamin localization during IFM development.

MotivationEnhancer-derived RNAs (eRNAs) and their fusion with protein-coding genes represent a crucial yet understudied layer of transcriptional regulation. eRNAs are typically expressed at low levels, which makes fusion events difficult to detect with conventional fusion detection tools. In addition, these tools are not designed to capture fusion transcripts arising from spatial proximity between distal regulatory elements and gene loci. Reads spanning such regions are also frequently filtered as mapping artifacts. As a result, computational approaches for systematically identifying spatially mediated enhancer-exon fusion transcripts remain lacking. MethodsWe developed ChiMER, a graph-based framework for detecting ChiMeric Enhancer RNAs from short-read RNA-seq data. ChiMER constructs splice graphs with chromatin contact information to introduce enhancer-exon edges and uses graph alignment to search for potential transcriptional paths. A ranking-based scoring module then prioritizes high-confidence events. Evaluations on simulated and real RNA-seq datasets show that ChiMER achieves higher sensitivity than conventional linear fusion detection methods while maintaining low false-positive rates. ResultsApplied to cancer cell line RNA-seq datasets, ChiMER identified multiple enhancer-exon chimeric transcripts, several associated with super-enhancer regions. Multi-omics analysis further shows that fusion transcripts occur in transcriptionally active regulatory environments and frequently coincide with strong R-loop signals, suggesting a potential role of RNA-DNA hybrid structures in facilitating long-range transcriptional joining events. Availabilityhttps://github.com/Candlelight-XYJ/ChiMER Contactyujia.xiang@outlook.com, xielinhai@ncpsb.org.cn

The MYC associated factor X (MAX) is the heterodimeric partner of the MYC paralogs (MYC, MYCN and MYCL). When deregulated, high level of the MYC paralogs contribute to all aspects of tumorigenesis and tumor growth. MAX can also heterodimerize with the MXD proteins, MNT and MGA. Heterodimerization and sequence specific DNA binding to the E-Box sequences at gene promoters is controlled by their heterodimerization with the MAX b-HLH-LZ. As a heterodimer with MAX, MYC proteins activate genes involved in cell metabolism, growth and proliferation whereas MXD proteins, MNT and MGA repress them. MAX can also bind to the E-Bos sequence as a homodimer. Being devoid of a transactivation domain it can act as an antagonist of the MYC/MAX heterodimers. Variants of MAX have been reported to be linked to cancer. These variants are either not expressed, inactivated or lead to missense mutations. This has led to the notion that MAX may have a tumor suppressor role. Here, we characterize three of those variants with missense mutations in the basic region, i.e. E32K, R35P and R35C. We analyzed their heterodimerization with the b-HLH-LZ of MYC and their DNA binding properties as homo-and heterodimers. The R35C variant b-HLH-LZ was found to have a markedly increased affinity for the b-HLH-LZ of MYC. We also observed that all three b-HLH-LZ variants have a lower affinity as homodimers for the E-Box than the WT. This was shown to lead to a preferential binding of all the heterodimeric b-LHLH-LZ to the E-Box. This effect is exacerbated in the case of the R35C variant. We argue that this preferential binding of MYC as heterodimers with these variants to E-Box sequences could contribute to tumorigenesis. Hence, our results suggest that, mechanistically, the MAX homodimer bound to the E-Box could act as a tumor suppressor. MATERIALS AND METHODSO_ST_ABSMolecular modelingC_ST_ABSThe open source version 1.7.6.0 of Pymol was used for modeling and molecular rendering [1]. The crystal structure of the MAX homodimer bound to the E-Box (1HLO [2]) was used as a template for the generation of the models. The variants were generated using the mutagenesis function in the wizard. The conformation of the K32 side chain was manually set in order to avoid introducing steric clashes with DNA. Protein expression and purificationThe cDNA, coding for the MAX b-HLH-LZ (Max* hereafter, residues 22-103, UniProt entry P61244-1) to which are added the GSGC residues in c-terminal, inserted in the pET3a vector was already available in the laboratory [3] and was used as a template to generate the plasmids with inserts coding for each of the mutants (E32K, R35C and R35P) through quick-change PCR with Q5 DNA polymerase and DpnI from New England Biolabs. The primers used were purchased from IDT DNA, their sequences are listed in Table S1. Sequence for each construct was confirmed by Sanger sequencing at the Plateforme de sequencage SANGER - Centre de recherche du CHU de Quebec - Universite Laval. The primary structure for the basic region of each construct is given in Fig. 2A. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=137 SRC="FIGDIR/small/715400v1_fig2.gif" ALT="Figure 2"> View larger version (41K): org.highwire.dtl.DTLVardef@1b05d5eorg.highwire.dtl.DTLVardef@1c1d692org.highwire.dtl.DTLVardef@ee469dorg.highwire.dtl.DTLVardef@15e0ba4_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 2.C_FLOATNO Structure schematics, specific and non-specific interactions dictating specificity and stability of binding of the basic region of MAX to the canonical (CACGTG) E-Box. A. Primary structure for the basic region of MAX and each of the variants. Positions making the most important contacts with the E-box are indicated by black arrows. Positions for the variants studied here are colored according to the Zappo colour scheme, following their physico-chemical properties: red for negative, blue for positive, magenta for proline and yellow for cysteine. B. The side chain (carboxylate) of E32 receives H-Bonds from the CA nucleobases in the leading strand (white carbon atoms). R35 and R36 make a salt bridges with phosphate groups while and the guanidino moiety of R36 makes a specific H-Bond with the nucleobase of the G in the strand of the reverse complement (cyan carbon atoms). C. The R35C mutation removes one non-specific salt-bridge at the interface of the complex. D. The aliphatic portion of the K side chain in the E32K variant is unable to accept the H-Bonds from the CA nucleobases and leads to the stabilisation of the complex and the helical structure of the basic region. E. In addition to removing a salt-bride, the Pro residue in the R35P kinks the path of the basic region, prevents the establishment of the specific H-Bonds mandatory for recognition of the E-Box and leads to unfolding of the helical state. C_FIG The MYC b-HLH-LZ (Myc*), the Max*WT b-HLH-LZ and its variants were expressed and purified as previously described [3,4] After lyophilisation, the b-HLH-LZs were kept at -20{degrees}C and solubilised in Myc buffer (50 mM NaCl, 50 mM NaH2PO4 pH 5.5) for Myc* or PBS for Max* at a final concentration of 1 mM before use. Circular dichroismAll circular dichroism (CD) measurements were performed on a Jasco J-810 spectropolarimeter equipped with a Peltier-type thermostat. The instrument was routinely calibrated using an aqueous solution of d-10-(+)-camphorsulfonic acid at 290.5 nm. Samples were prepared as follows: Max* (either WT or a variant) was diluted in 100 {micro}l 2X CD buffer (40 mM KCl, 11.4 mM K2HPO4, 28.6 mM KH2PO4, pH 6.8) and the volume adjusted to 106 {micro}l with PBS. 10 {micro}l TCEP 16 mM were added, and the volume further adjusted to 192 {micro}l with ddH2O before samples were incubated overnight at room temperature. After reduction, Myc* was added and the volume adjusted to 198 {micro}l with Myc buffer (Na2HPO4 0.95 mM, NaH2PO4 49.05 mM, 50 mM NaCl, pH 5.5). The DNA complexes were prepared as follows. After a 10 minutes incubation of the protein samples at room temperature, 0, 1 or 2 {micro}l of 2 mM of specific or non-specific DNA duplexes in 10 mM Tris pH 8.0 were added and the volume adjusted to 200 {micro}l with 10 mM Tris pH 8.0. The strands of the specific probe were: 5-ATT ACC CAC GTG TCC T*AC-3 and 5-GTA GGA CAC GTG GGT* AAT-3 (with the E-box sequence underlined) and the non-specific probe: 5-ATT ACC TCC GGA TCC T*AC-3 and 5-GTA GGA TCC GGA GGT* AAT-3 (Integrated DNA Technologies). Samples were further incubated for 10 minutes at room temperature and transferred to a 1 mm path length quartz cuvette. All spectra were recorded from 250 to 195 nm at 0.1 nm intervals by accumulating 10 spectra at 25 {degrees}C. Thermal denaturations were recorded at 222 nm from 5 to 95 {degrees}C at a heating rate of 1 {degrees}C/min. CD signal for spectra and thermal denaturations was corrected by substracting the signal from corresponding spectra or thermal denaturation either for buffer alone or the appropriate DNA duplex. CD signal was then converted to mean residue ellipticity using the following formula [5]: [{theta}] = {delta} {middle dot} MRW/(10{middle dot}c l) where [{theta}] is the mean residue ellipticity in deg {middle dot} cm2 dmol-1, {delta} is the CD signal in millidegrees, MRW is the mean residue weight, c is the concentration in mg/ml and l is the pathlength in mm. For the heterodimers, the concentration used was the sum of Max* and Myc* and the MRW was determined using a weighted average.

This retrospective study aims to explore the interactive effects of biological maturation and relative age effect (RAE) on talent identification. 56 male elite soccer players matched for chronological age (15.08{+/-}0.41 years) were studied. Test items included anthropometry (height, body mass, sitting height, leg length, BMI and Quetelet index), physiology (power, speed, agility, speed endurance and aerobic performance), soccer-specific skills (passing, shooting and dribbling), psychology (achievement motivation, orientation and resilience) and biological maturation (APHV) tests. The test results were analyzed independent sample t-test, Pearson correlation analysis, and stratified regression. Conclusion: Biological maturation significantly influences anthropometry (height, weight and Quetelet index), lower limb explosive, and speed (single-leg jump, standing triple jump, and 30-m sprint) in U16 male elite soccer players in Shanghai. The relative age effect shows no significant impact on talent selection indicators, which is attributed to the accumulated training load effect. The mechanisms of biological maturation and RAE in youth soccer talent selection are distinct and operate independently.